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Abstract 


The transport and mixing characteristics of a large sample of air parcels within 
a mature and vertically sheared hurricane vortex is examined. Data from a high- 
resolution (2 km grid spacing) numerical simulation of “real-case” Hurricane Bonnie 
(1998) is used to calculate Lagrangian trajectories of air parcels in various subdomains 
of the hurricane (namely, the eye, eyewall, and near-environment) to study the degree 
of interaction (transport and mixing) between these subdomains. It is found that 1) 
there is transport and mixing from the low-level eye to the eyewall that carries high- 
6 e air which can enhance the efficiency of the hurricane heat engine; 2) a portion of 
the low-level inflow of the hurricane bypasses the eyewall to enter the eye, that both 
replaces the mass of the low-level eye and lingers for a sufficient time (order 1 hour) to 
acquire enhanced entropy characteristics through interaction with the ocean beneath 
the eye; 3) air in the mid- to upper-level eye is exchanged with the eyewall such that 
more than half the air of the eye is exchanged in five hours in this case of a sheared 
hurricane; and 4) that one-fifth of the mass in the eyewall at a height of 5 km has an 
origin in the mid- to upper-level environment where 6 e is much less than in the eyewall, 
which ventilates the ensemble average eyewall 0 e by about 1 K. Implications of these 
findings to the problem of hurricane intensity forecasting are discussed. 



1 Introduction 


Vortices exhibit the well known property that mixing between their cores and near envi- 
ronment is strongly suppressed. Examples are manifold, including vortex rings (Vladimirov 
and Tarasov 1979), mesoscale ocean vortices (Robinson 1982), emergent vortices in quasi- 
geostrophic and two-dimensional turbulence (Carnevale et al. 1991; McWilliams et al. 
1994), the wintertime stratospheric vortex (McIntyre 1989), and tropical cyclones (TCs) 
(Willoughby 1998). The details of these flows are quite different, yet the tendency for 
vortices to act as mixing barriers to environmental fluid is believed central to both their 
emergence and long lifetimes in complex flows (e.g., McWilliams 1984; Mizuta and Yoden 
2001). This vortex robustness can be traced to the quasi-elastic behavior of vortices, which 
results from having mechanisms that can oppose disturbances that are either axisymmetric 
[via centrifugal restoration (Howard and Gupta 1962; Charney 1973)] or nonaxisymmetric 
[via realness of discrete spectra (Vladimirov and Tarasov 1980) or via a tendency to axisym- 
metrize (Melander et al. 1987; Montgomery and Kallenbach 1997; Bassom and Gilbert 1999; 
Reasor et al. 2004; Schecter and Montgomery 2004)] . If the robustness of a vortex core were 
to be reduced, however, this could have important consequences on the maintenance of that 
vortex. 

In TCs, the manifestation of some compromise in the robustness of the vortex is the 
stirring 1 of air (and attendant material properties such as angular momentum or equivalent 
potential temperature) between the various domains of the storm, i.e. the eye, the eyewall, 
and the surrounding environment. Among the key atmospheric mechanisms that have been 
identified to produce mixing between TC domains are transport by frictionally-induced in- 
flow that slips under the eyewall (Persing and Montgomery 2003, hereafter PM03; Braun 

1 We follow the definitions found in Haynes (2002). “Transport” is to carry air with various properties 
from their source to their sinks. “Stirring” is to bring air with different properties into closer proximity with 
each other. Stirring which proceeds to the point that the originating characteristics of the air can no longer 
be distinguished is called “mixing”. “Mixing” will also be used in general discussion where the distinction 
between these processes is not important. 
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2002) ; stirring by barotropic /baroclinic (inflectional) instabilities that break down the high 
potential vorticity annulus comprising the eyewall (Schubert et al. 1999; Kossin and Schu- 
bert 2001; Nolan and Montgomery 2002); and stirring by an imposed vertical wind shear 
(Frank and Ritchie 2001; Zehr 2003; Emanuel et al. 2004). Each of these processes has been 
shown capable of significantly altering the intensity of TCs (Wu and Braun 2004; Kossin 
and Eastin 2002; Bender 1997; Wang 2002; Emanuel 1989, 1997; Montgomery et al. 2005; 
Prank and Ritchie 2001; Wong and Chan 2004). While signifying a compromise in the ro- 
bustness of a vortex, these mixing processes can alter the thermodynamics of the hurricane 
eyewall, whether by enhancement or dilution of relevant thermodynamic properties. Recall 
for example, that the uptake of moist enthalpy from the underlying ocean and the ascent 
in the eyewall of the TC represent the isothermal and moist adiabatic legs of an idealized 
Carnot cycle, respectively. An enhanced intake of energy to the Carnot cycle at this stage 
suggests an increase in useful work (enhancement of the kinetic energy field) and ultimately 
an increase in the loss of energy to frictional dissipation at the ocean surface (Emanuel 1986, 
1995; PM03) . 

Evidence in support of the hypothesis that instabilities in the TC eyewall, and their finite 
amplitude coherent structures (such as illustrated in Schubert et al. 1999 and Montgomery 
et al. 2002), contribute positively to the intensity of the storm was presented in Emanuel 
(1989, 1997) and PM03. In Emanuel (1997), the breakdown of the vortex sheet comprising 
the eyewall of a developing storm was shown to be essential for achieving the maximum 
potential intensity predicted by the Carnot theory (hereafter E-MPI). In PM03, eyewall 
mesovortices bordering the low-level eye and eyewall region were spawned by a (Kelvin- 
Helmholtz) shear instability on the annulus of high toroidal vorticity associated with the 
secondary circulation of the hurricane. A characteristic of these coherent structures, not 
anticipated by Emanuel (1989, 1997), was their ability to persistently flux high entropy 
air from the low-level eye to the eyewall. The temporal eddy and mean transport of high 
entropy air from the low-level eye to the eyewall was shown to provide supplementary “fuel” 
to the heat engine, in addition to what is drawn from the ocean underneath and outside the 
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eyewall. By utilizing this extra fuel, the hurricane can obtain a “superintense” state (PM03; 
Montgomery et al. 2005), in which the storm exceeds E-MPI. Even if a storm in nature 
is not superintense, PM03 suggested that the superintensity mechanism may still operate 
to partially mitigate the negative influences associated with upwelling cooler water, cold 
downdrafts, environmental vertical shear, or the interaction with mid-tropospheric dry air. 

A complementary situation to superintensity is the problem of transport and stirring 
of mid-tropospheric dry, low-entropy air from outside the rain area into the eyewall region. 
This is a form of entrainment that occurs on the vortex scale. Previous work (Simpson and 
Riehl 1958; Riehl and Malkus 1961; Gray 1968; Emanuel et al. 2004) has suggested that 
such an exchange process, which can result from an imposed vertical wind shear, can dilute 
the heat content of the eyewall thereby weakening the intensity of the TC. This process was 
coined “ventilation” by Simpson and Riehl (1958). In practical situations, ventilation may 
occur as vertical wind shear impinges on a storm or, alternatively, as a storm approaches 
a middle- and upper-tropospheric air mass with low entropy. To give sufficient warning 
to coastal communities, forecasters must decide if the TC core will be adversely affected 
by this interaction 2 . In general, ventilation requires dry air to breach the intense vorticity 
of the eyewall. Idealized modeling studies suggest, however, that this can only occur via 
nonlinear breakdown of the annulus of high potential vorticity comprising the eyewall and 
the attendant stirring between the environment and eyewall air masses (e.g., Schubert et al. 
1999). Ventilation is unlike the superintensity mechanism that provides “extra” fuel from 
the low-level eye region. Processes that entrain dry air into the TC core will be referred to 
here as the “anti-fuel” problem. 

Little is known quantitatively about the efficiency of the fuel or anti-fuel mechanisms 
in real-case TCs. In this paper we take the first step in a more systematic study to test 
and quantify these ideas. The real-case numerical simulation of Hurricane Bonnie (1998) by 
Braun et al. (2005) employing the Pennsylvania State University-National Center for At- 
mospheric Research (PSU-NCAR) non-hydrostatic fifth generation mesoscale model (MM5 

2 In this paper, the “core” of the TC refers to the eye and the high- vorticity region of the eyewall. 
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V3.4; Dudhia 1993; Grell et al. 1995) provides a convenient starting point for our examina- 
tion. By using a large sample of air parcels seeded throughout the storm we build upon the 
study of Braun (2002) and compute several thousand three-dimensional Lagrangian trajecto- 
ries for air parcels seeded throughout the eye, eyewall, and surrounding regions of Hurricane 
Bonnie. The large number of Lagrangian trajectories is used to gather a “census” of the 
behavior and thermodynamic properties of air parcels in various regions of the simulated 
hurricane. 

Although Bonnie represents only a single case, and results gleaned from it cannot easily be 
generalized, we believe the methodology developed herein is nonetheless useful for laying the 
groundwork for a more systematic study of the hurricane mixing phenomena using idealized 
TC flow configurations. Our objective here is threefold: Identify the nature of the mixing 
processes that occur within this vertically-sheared storm; quantify the efficiency of the mixing 
between the eye, eyewall and enviroment; and relate the current findings on mixing within the 
vortex to the heat transport that maintains the vortex at approximately constant intensity 
during a period of near-constant vertical shear. 

The outline of the paper is as follows. Section 2 describes the model setup and provides an 
overview of the Hurricane Bonnie simulation. Section 3 describes the Lagrangian trajectory 
seeding strategy and provides results illustrating the mixing that occurs between various 
regions of the model hurricane. Figure 1 provides a simplified view of these mixing processes 
in the context of a hurricane that experiences a moderate degree of vertical wind shear 
(10 m s -1 /10 km). Section 3.1 focuses on the low-level eye/eyewall interaction, including 
evidence of the superintensity mechanism (trajectory class I in Fig. 1); Section 3.2 focuses on 
trajectories in the boundary layer inflow layer and a quantification of the low-level eye mass 
replacement (class II in Fig. 1); Section 3.3 focuses on the upper-level eye/eyewall interaction 
and investigates the eye “containment” hypothesis proposed by Willoughby (1998) (class III 
in Fig. 1); and Section 3.4 focuses on the interaction of the midlevel environment with the 
storm and how the environmental vertical shear may serve to ventilate the eyewall and/or 
eyewall region (class IV in Fig. 1). Section 4 concludes with a summary of the results and 
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objectives for future research. 


2 Model setup and overview of the Bonnie simulation 

Braun et al. (2005) presented an in-depth analysis of a numerical simulation of Hurricane 
Bonnie using the PSU-NCAR MM5 (version 3.4). We use the same simulation in this study 
but focus on the mixing characteristics that can be deduced from the use of a large number 
of Lagrangian trajectories. 

2.1 Experimental design 

Summarizing the numerical setup, four nested grid domains of increasing horizontal resolu- 
tion were used. Due to computational limits, the simulation was carried out in two steps. 
First, the two coarser meshes were simulated, using a 36-km horizontal grid spacing and 
91 x 97 x 27 grid points in the x , y, and z directions for the first grid, and 160 x 160 x 27 
grid points at 12-km horizontal grid spacing on the second grid. The simulation on these 
two outer grid meshes was started at 1200 UTC 22 August 1998 and run for 36 hours, with 
model output saved every hour. A modified version of the Blackadar planetary boundary 
layer scheme is employed in which surface roughness is calculated for momentum, tempera- 
ture, and moisture following Garratt (1992) and Pagowski and Moore (2001). Cloud effects 
are parameterized on the coarse mesh in two ways: the Grell cumulus scheme (Grell et al. 
1995) is used to represent unresolved convective processes while the Goddard microphysics 
scheme (Tao and Simpson 1993; McCumber et al. 1991) is used to represent grid-scale cloud 
microphysical processes. 

A higher-resolution simulation was then carried out using the 1-h output from the 36- 
km and 12-ktn grid results to provide initial and boundary conditions. The model grids 
consisted of a 6-km grid with 225 x 225 x 27 grid points and a 2-km nest with 226 x 226 x 27 
grid points. The simulation on these high-resolution grids was started 6 h into the coarse- 
resolution forecast to allow for model spin up on the 12-km grid, and was carried out for 
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30 h until 0000 UTC 24 August. The 2-km grid was moved periodically to keep it centered 
over the storm. Model physics were similar to the coarse-grid runs except that no cumulus 
scheme was used. Further details of the model setup, initialization, and model physics are 
provided in Braun et al. (2005). For this paper, t = Oh will be considered the start of the 
fine mesh simulation, i.e. 1800 UTC 22 August 1998. 

This study focuses on a 5 h period from t = 15 — 20 h, 0900-1400 UTC 23 August 1998. 
The simulated storm during this time undergoes an intensifying trend (Fig. 2). The minimum 
surface pressure p sfc decreases steadily from 961 mb at 14 h to 951 mb at 20 h. Maxima of 
tangential wind v sfc at 0 — 40 m during this time vary between 50 and 55 m s -1 , although 
there is a slight increasing trend in the 6 h period. The maximum azimuthally averaged 
tangential wind v Bfc increases slightly from approximately 40 to 44 m s -1 . Comparison with 
NHC Best Track for Hurricane Bonnie at (0600; 1200; 1800) UTC shows p sic intensifying 
(960; 958; 955) mb and wind speeds increasing (49; 51; 51) m s -1 . 

2.2 Summary of simulated structure 

A complete summary of the simulated structure of the Bonnie simulation can be found in 
Braun et al. (2005). Several features of particular interest to the current study, however, 
are summarized below. 

The system-scale structure of the storm is primarily influenced by northwesterly envi- 
ronmental deep layer (850-200 mb) vertical shear. Because vertical wind shear has been 
identified (DeMaria and Kaplan 1994) as a leading predictor of storm intensity, we should 
anticipate that shear will be the most important factor in determining simulated intensity. 
Simulated radar reflectivity (shown in Fig. 6 in Braun et al. 2005) exhibits an asymmet- 
ric storm structure consistent with the influence of northwesterly vertical wind shear. The 
highest reflectivities appear southeast, east, and north of the storm center, which is down- 
shear and left of the 850-200 mb vertical shear vector calculated about the storm center. 
The vertical wind shear has a magnitude of approximately 12 m s -1 between z = 1 and 
^ = 12 km; observations from other TCs have shown that at this magnitude, vertical shear 
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effects dominate subsequent storm evolution (e.g., Black et al. 2002). The boundary layer 
is characterized by inflow (Fig. 3a) in all quadrants outside the radius of maximum wind 
(RMW), but is stronger southeast, east, and north of the vortex center. 3 Generally, the peak 
near-surface inflow is confined below z = 1 km, decreasing to 2 = 750 m near the RMW. In 
all regions of the low-level eye except north of the center, radial outflow 'prevails and reaches 
in excess of 5 m s -1 . Radial flow into and out of the core of the storm can be found at other 
levels in association with the environmental shear profile. 

Just above the surface inflow layer (z == 1.1 km), the region outside the RMW exhibits 
radial outflow to the west and inflow to the east of the storm center. Inside the RMW, radial 
outflow exists directly downshear and left of the deep layer mean vertical shear vector. The 
outflow inside the RMW at these two levels (z = 453 m and 1.1 km) is associated with 
the convective updrafts which preferentially form in the downshear side of the eyewall, as 
described in Braun et al. (2005). 

In the middle and upper troposphere (z = 5.6 and 9.9 km) radial outflow occurs south, 
east, and northeast of the storm center. Inside the RMW, the outflow is particularly strong 
(u > 10 m s -1 ) directly south of the storm center. This region was shown in Braun et al. 
(2005) to be a favorable area for the initiation of eyewall convective “hot towers” . 

The azimuthal mean storm structure (Fig. 4) exhibits the classical outward sloping eye- 
wall of a mature hurricane with maximum mean tangential wind v of 55 m s -1 and mean 
vertical velocity w of 1.2 m s -1 . The RMW near the surface at this time (t = 15 h) is 44 
km and gradually moves inward with time to 38 km at t = 20 h (not shown). Maximum w 
increases from 1.2 to nearly 2.3 m s -1 through the period t = 15 — 20 h. 

The distribution of 9 e air (Fig. 4c, following Bolton 1980) suggests a reservoir of high 
9 e in the boundary layer of the eye that is characteristic of mature hurricanes (Willoughby 
1998; Zhang et al. 2002). In the eye, the maximum in 0 e is greater than 376 K near the 

surface and reaches a minimum at z = 4 km. 

3 The method for computing the storm center is described in the appendix. 
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3 Use of trajectories to study transport, stirring, and 


mixing processes within the vortex 

Each of the following sections will illustrate the interaction between various regions of the 
model hurricane using distinct sampling strategies for the initial locations of trajectories. 
The strategies are detailed separately in each section. A simplified view of a moderately 
sheared hurricane is shown in Fig. 1 and also describes the transport processes that we 
are interested in analyzing in this study. The calculation of trajectories is performed by 
post-processing the model output data, which has a temporal output increment of three 
minutes. Trajectory positions are computed every 10 seconds by the technique described in 
the Appendix A. 

3.1 The low-level eye/eyewall interaction 

The trajectory sample 

For this section, we examine the interaction between the low level eye and the eyewall using 
forward trajectories (class I from Fig. 1) seeded in the eye below the eye inversion. The 
trajectories are seeded at 2 == 453 m in the eye at t — 900 min = 15 h from the start of the 
inner grid simulation and are calculated forward for 5 h from this time. The starting locations 
are on a rectilinear grid with a 2 km spacing in both x and y, matching the specification of 
the model grid. 

The eye is distinguished from the eyewall using the following axisymmetric criteria. The 
eyewall is defined here as the region in the vicinity of the RMW with azimuthally averaged 
vertical velocity w > 0.25 m s _1 and axisymmetric mean relative humidity 7 i > 90% . The 
low-level eye (below z = 1 km) is defined as the region inside the RMW and near the storm 
center with 9 e > 370 K. The eye at z = 1.1 km is defined as the region inside the RMW and 
near the storm center with w < 0.25 m s -1 . The eye at mid- to upper-levels (z = 5.5 km and 
above) is defined as the region inside the RMW and near the storm center with w < 0.25 
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m s“ 1 and H < 90% . These definitions using axisyrnmetric properties help distinguish the 
eyewall from short-lived cumuli in the eye itself. Although it would be preferable to apply 
the same definition for the eye at z = 1.1 km as it is at mid- to upper-levels, the relatively 
higher values of 7i at 1.1 km (~ 90-92% ) does not provide a clear distinction between the 
mean eye and eyewall and thus precludes the use of a humidity discriminant at this level. 
The shear-induced asymmetry in the storm lowers the axisyrnmetric relative humidity below 
saturation, hence the value of H = 90% is chosen as the humidity discriminant for defining 
the mean eyewall (cf. Fig. 4). 

Transport from the low-level eye to the eyewall 

A small sample of the trajectories initiated in the eye that enter the intense eyewall updraft 
are shown in Fig. 5. After first encountering the eyewall, these seven trajectories encounter a 
downdraft at z ~ 1 — 2 km and approach the center of the storm before becoming entrained 
in the eyewall updraft. At upper levels (z > 7 km), six of these trajectories detrain into the 
eye. Vertical velocities of the eyewall updraft encountered by these trajectories range from 
1 — 2 m s" 1 below z = 3 km to as much as 10 m s -1 above this height (not shown). The 
initial 9 e associated with the trajectories is between 372 and 376 K, confirmation that these 
trajectories are originating from the reservoir of relatively high 9 e located in the low level 
eye (9 e greater than 370 K extend up to 2 = 750 m at the storm center; cf. Fig. 4c). As 
the trajectories rise in the eyewall updraft, the associated 9 e decreases approximately 5 K. 
The vertical profile of 9 e along the trajectories becomes relatively constant with height above 
z = 2 km. These parcels of air, originating in the low-level eye, retain a degree of warmth 
relative to other eyewall parcels, as is seen by plotting asymmetric Q e along the trajectories 
(Fig. 6a). Asymmetries in this paper are defined by subtracting the axisyrnmetric mean 9 e 
field from the complete 9 e field. Below z — 3 km, asymmetric w (Fig. 6b) is both negative 
and positive, the former being associated with the downward and inward motion exhibited 
in Fig. 5b. Asymmetric w is positive between z = 3 km and 7 km, with values as high as 
10 ms -1 . The asymmetric pictures of 9 e and w suggest that air stirred out into the eyewall 
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from the low-level eye retains some degree of relative buoyancy to other “default” parcels of 
eyewall air. This relative buoyancy can be associated with the “hot towers” identified by 
Braun (2002) and Smith et al. (2005). 

Consider now the complete subset of trajectories that make their way out to the axisym- 
metric eyewall (Fig. 7). There is an overall tendency for these trajectories to be stirred out 
into the eyewall downshear and left of the shear vector, and stronger tendency is shown for 
those trajectories (dark dots) that are stirred out near the base of the eyewall. Of the 551 
trajectories which are seeded inside the eye at this level, 313 trajectories, or 56.8% , are 
mixed out to the eyewall within 5 h. Although the “escape” points shown are located all 
around the RMW azimuth ring, the majority of trajectories (Fig. 7b) encounter the eyewall 
at locations southeast, east, and northeast of the storm center. These findings are consistent 
with the results in Braun et al. (2005), which show that the main eyewall updrafts occur 
downshear and left of the mean layer vertical shear vector. 

Evidence of “superintensity” mechanism 

Persing and Montgomery (2003; PM03) defined superintensity as a state where hurricane 
intensity exceeds the predicted maximum intensity of Emanuel’s (1995) maximum potential 
intensity (E-MPI) theory. The mechanism identified by PM03 in their axisymmetric mod- 
eling was transport of air from the low-level eye to the eyewall. The low-level eye air was 
enhanced in 0 e \ thus, its introduction to the base of the eyewall represented an additional 
source of heat for the hurricane heat engine cycle. In E-MPI theory, the only source of heat 
for driving the hurricane circulation is the ocean beneath the eyewall through a wind-induced 
heat exchange. 

The representative sample of trajectories shown in Fig. 5 provides evidence of the super in- 
tensity mechanism operating in the MM5 Bonnie simulation. Once these sample trajectories 
are stirred out to the eyewall, they are lifted by the eyewall updrafts. The form of the inter- 
action between the high-# e trajectory particles and the eyewall air is an incomplete reduction 
of trajectory 6 e to the mean eyewall value, retaining a moderate localized enhancement of 
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the eyewall updraft associated with the trajectory (Pig. 6). Since 0 e is approximately con- 
served, 4 a reduction of 6 e following the trajectory upon introduction to the eyewall must 
be associated with an increase in 9 e in the surrounding air, thus warming the axisymmetric 
mean 6 e by mixing. 

The stirring of high-# e air from the low-level eye into the eyewall is corroborated by 
instantaneous r — z cross sections of 0 e along radial legs in the southeast quadrant of the 
storm (not shown). As shown in Braun et al. (2005), the eyewall updrafts originate in this 
quadrant of the storm, which corresponds to the downtilt-right side of the eyewall in relation 
to the environmental vertical shear vector. These cross sections of 9 e exhibit a tongue of 
high-0 e air in the low level eye extending into the base of the eyewall updrafts [see an example 
of this structure in Fig. 12 of Braun (2002)]. This suggests that this warmer 0 e air is being 
introduced to the eyewall updrafts as they form. 

For the case of the Bonnie simulation, the sea surface temperature is 303 K, the latitude 
is 24.2° N, and the ambient surface pressure is 1010 mb. Since there is no dissipative heating 
in the numerical simulation, Emanuel (1995) represents the appropriate derivation of E-MPI 
for comparison. A reasonable outflow temperature following the technique of PM03 of the 
axisymmetric circulation is 210 K. To supply the values of surface relative humdity, we choose 
RH=80% as an expedient. An analysis of Ck/Co for the surface flux parameterization for 
this model run exhibits a spatial variability for the ratio, ranging from 0.65 for the average 
over the complete inner domain to 0.35 at the radius of maximum winds. With these values, 
E-MPI predicts a value of maximum v at the top of the boundary layer of 52.2 m s -1 for 
C k /C D = 0.65 to 38.2 m s -1 for Ck/Co = 0.35. The strongest axisymmetric tangential 
winds of 55 m s -1 found at (r = 48, z — 0.8 km) in the Bonnie simulation are near to or 
much greater than E-MPI. The simulation presented here is thus plausibly superintense. 

4 The MM5 model is not designed to conserve 6 e of Bolton (1980). Both 9 e and MM5, however, are 
designed in reference to reality, so approximate conservation can be anticipated (and is rougly observed 
in the post analyses) in the absence of ice processes, surface interaction, or strong sub-grid scale diffusive 
mixing. 
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3.2 Eye inflow trajectories and replenishment of the 9 e reservoir 

The trajectory sample 

If the eyewall updrafts are consuming some of the high 9 e air within the low-level eye, then 
there must be some mechanism for replenishing that high 9 e reservoir. To examine this 
replenishment process, in this section, we investigate forward trajectories (class II of Fig. 1) 
seeded in the inflow layer converging on the storm center. The question to be addressed at 
this point is: what percentage of trajectories are transported underneath the eyewall updrafts 
and into the hurricane eye, as opposed to following the standard ascent path through the 
eyewall updrafts? 

Seed points are located horizontally on a radius vs. azimuth grid. Radial seed spacing is 5 
km, between 50 < r < 180 km. Twenty-five azimuthal seed points exist at each seed radius. 
The vertical seed levels used are z = 40, 121, and 244 m. This seeding strategy results in 625 
trajectories calculated at each vertical seed level. The trajectories are calculated forward 3 
h from a seed time of t = 15 h into the 2 km simulation. 

To ensure that trajectories seeded at or near the lowest model level (z = 40 m) do not 
cross the lower grid boundary, it was necessary for us to extrapolate the velocity data to 
a near-surface height; we chose z = 10 m. Here we follow the logarithmic boundary layer 
as observed by Powell et al. (2003) 5 and assign the magnitudes of the u- and u- velocity 
components to be 85% of their values at z = 40 m. We assign the vertical velocity w to 
vanish at z = 10 m in order to ensure that no trajectories descend below 10 m. , 

Low-level eye mass replacement 

For a trajectory of class II to be transported into the eye, we required the 9 e value associated 
with a particular trajectory to exceed 370 K below a vertical height of z — 2 km. This 
threshold value is chosen based on the values of 9 e in the high-# e reservoir of the low level eye 
(Fig. 4c) . By this criteria, a significant fraction of class II trajectories (i.e., the boundary layer 
5 We assume that the parameterization of the surface layer in the model follows similar assumptions. 
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inflow) spends at least some time in the eye (Fig. 8). Considering trajectories originating 
just outside the eyewall (r < 100 km), close to half of low level (z — 40 m seeding; Fig. 8c) 
air is found to slip under the eyewall and enter the eye, with probability decreasing with 
height to about 20% at z = 121 m. At further radii (r > 100 km), turbulent redistribution 
between levels within the boundary layer blurs the distinction between'these levels. 

We now seek trajectories that can serve as prototypes for the two primary consequences 
for class II trajectories that reach the eyewall: either to enter the eye for some period of 
time or to be immediately drawn into the eyewall updraft without encountering the eye. 
The approach for constructing these prototypes will eliminate a large number of borderline 
cases, since these trajectories typify neither consequence. For the “eye inflow” prototypes, 
consider those trajectories of class II defined above which exceed 0 e = 372 K at some point 
in the calculation. For the “eyewall ascent” prototypes, consider those trajectories of class 
II with maximum 6 e between 369 K < 0 e < 370 K (Fig. 9) . The amount of time eye-inflow 
prototype trajectories spend in the eye can vary greatly, where this time is defined as the 
period spent with 0 e > 370 K just prior to entrainment into the eyewall updraft. The amount 
of time required for a trajectory to acquire a substantial increase in 6 e (3 K) can be as short 
as 15 min (Fig. 10), although 40-60 min is typical. The behavoir described here is sufficient 
to explain how the eye entropy reservoir is maintained; the mass to replace air entrained into 
the eyewall can be adequately accounted for while it is demonstrated that the 6 e of these 
parcels can be quickly enhanced through exchange with the underlying ocean. 

Superintensity mechanism maintained 

The “eye inflow” trajectories (Fig. 9a, c) show clear evidence of the superintensity mechanism 
occuring. All trajectories converge toward the center near the surface (z < 200 m) and 
approach to within 15-25 km of the storm center as they are stirred into the eye. When they 
are stirred out to the eyewall and lifted by the eyewall updrafts (Fig. 9c), the interaction 
between the high-f? e trajectory particles and the eyewall air is apparent. The 0 e of the 
trajectories is not maintained at the values found in the low-level eye, but decrease partially 
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to the value of the surrounding eyewall air as a result of diffusive mixing. Conversely, the 
“eyewall ascent” trajectories (Fig. 9b, d) show a slight warming initially in 9 e as they ascend 
the eyewall updraft. 

After exiting the eye, most of the “eye inflow” trajectories exhibit 6 e decreases of as much 
as 4 K (Fig. 11), as a result of mixing with relatively cooler eyewall air. The “eyewall ascent” 
trajectories, in turn, exhibit increases in 6 e - up to 3 K — resulting from their interaction 
with the higher 9 e air coming from the low-level eye. 

As the moist entropy of the low level eye is stirred out to the eyewall, the reservoir of high 
9 e air in the low level eye must be replenished. This is achieved primarily via the convergence 
of air toward the storm center in the low-level inflow layer as shown above (Fig. 9a). 

3.3 The upper-level eye/ eyewall interaction 

The trajectory sample 

To study the degree of isolation of the eye from the eyewall in the mid- to upper troposphere, 
we examine in this section forward trajectories' (class III of Fig. 1) seeded within the eye 
above the inversion. The procedure follows section 3.1 except that seeding at heights of 1.1, 
5.6, and 9.9 km will also be considered. The latter two are taken to represent' air parcels 
from the mid- to upper level eye, while 1.1 km is more representative of the lower level eye. 

The “containment vessel” hypothesis reviewed 

For trajectories originating inside the eye in the middle and upper troposphere, we are 
interested in whether the eye behaves as a “containment vessel” . Proposed by Willoughby 
(1998; W98), the containment vessel hypothesis posits that the air contained inside the eye 
above the eye inversion in intense tropical cyclones has remained inside the eye since the 
eye formed. Using flight-level and dropsonde data, W98 identified eye soundings typical 
of intensifying tropical cyclones. The main features of these soundings included an eye 
inversion between 900 and 850 mb and dew point depressions of approximately 10 K above 
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the inversion. The rise and/or fall of the inversion results from a balance between low level 
moist inflow, which is brought into the eye to compensate for the moist low level eye air 
lost to the eyewall updrafts (as shown previously in section 3.2), and warm, dry subsidence 
above the inversion induced by the the loss of low-level mass to the eyewall updrafts. W98 
describes the convective eyewall updrafts as “heat pumps” which do work on the eye by 
drawing moist air out of the low level eye and hence force thermally indirect descent above 
the eye inversion. 

The eye soundings typical of weakening hurricanes (W98) showed eyewall moisture being 
mixed into the eye at mid- to upper-levels. In contrast to intensifying storms, the weakening 
case is characterized by a more gradual, “Y-shaped” inversion. The base of the inversion in 
the weakening cases is typically found (W98) between 850 and 700 mb. 


Evidence of upper-level eye mass recycling 

The MM5 Bonnie simulation described in Braun et al. (2005) and herein appears to be 
more typical of the soundings in W98 in which eyewall moisture is being stirred into the eye 
(Fig. 12). The center eye sounding follows a saturated adiabat between 900 and 650 mb, 
with the exception of a shallow isothermal layer between 850 and 800 mb. The air between 
850 mb and the ocean is moister, with a fairly constant dew point depression of 1 K. Above 
the isothermal layer,.dew point depressions range between 2 and 6 K. 

Trajectories seeded in the eye at the z = 5.6 and 9.8 km levels suggest that more than 
half (59% at 5.6 km, 76% at 9.8 km) of the eye air is removed from the axisymmetric eye to 
the axisymmetric eyewall (using the same criteria as above) over the 5 h period of study. At 
z = 5.6 km, the exit locations for trajectories (Fig. 13) tend toward the south east sector of 
the eye, consistent with the radial outflow shown at these azimuths (Fig. 3c). At z = 9.8 km, 
the pattern is the same, but shifted toward southern azimuths (Figs. 14, 3d). Braun et al. 
(2005) has associated this radial wind signature with the outflow from strong updraft towers 
which form on the southern side of the eyewall. 

Figure 15 shows the percentage of eye trajectories seeded at low- to mid-levels (z = 453 m, 


15 



1.1 km, and 5.6 km) which are stirred into the eyewall. In general, the percentage increases 
with increasing seed radius, and the percentage of trajectories seeded at lower levels (z = 
453m and 1.1 km) that are being stirred into the eyewall is approximately 10-20% higher 
than what is shown at the 2 = 5.6 km seed level. A large fraction (45% or more) of air 
from lower levels (z = 453 m and 1.1 km) of the eye encounter the eyewall within 5 h to 
within r = 15 km of the center, evidence for a large degree of mass replacement of the eye. 
Although we have not comprehensively documented the mid- to upper-level eye replacement 
(which would require a longer time interval), the process is illustrated in part by the sample 
of seven class I trajectories (Fig. 5), where air rising through the eyewall detrain to the eye 
above z = 6 km. 

3.4 Ventilation of mid-level environment to the TC 

3.4.1 The trajectory sample 

To examine the degree to which dry, mid-level environmental air intrudes upon the eyewall 
(ventilation), for this section we investigate backward trajectories seeded across the eyewall 
and determine the sources of that air. All four classes of trajectories found in Fig. 1 can 
in principle be found in such a sample, and the objective will be to distinguish class IV 
trajectories from the sample and examine their thermodynamic characteristics. 

The seed locations are at z = 5 km with equal radial/azimuthal spacing between 36 < 
r < 50 km. There are 50 seed points at each radius and the radial spacing is 2 km, for a 
total of 400 trajectory seeds. The seed time used is t = 20 h and the backward calculation 
is carried out for 5 h. Throughout the seeding annulus, w > 0.25 m s _1 . 

The process for distinguishing different classes of trajectories is as follows. For each back- 
trajectory, W > 0.25 m s -1 at the seed (final) time, by construction. Each trajectory may 
be followed backward in time from the seed location to determine the first location where w 
is no longer greater than > 0.25 m s -1 , which we define as the eyewall entry location. It is 
possible that one or more trajectories might remain in the axisymmetric eyewall throughout 
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the trajectory calculation, but it was found that all trajectories had origins outside the 
eyewall through the 5 h calculation period. Trajectories with eyewall entry points below 
z — 1.5 km are considered in either class I or II, and no effort is made in the rest of this 
section to distinguish class I from class II. Trajectories with eyewall entry above z — 1.5 km 
but radius interior to the radius of maximum updraft are considered class III. For the rest 
of the trajectories (eyewall entry above z = 1.5 km and outside the radius of maximum 
updraft), two modes were identified (Fig. 16a). The first mode of trajectories were involved 
in updrafts just exterior to the axisymmetric eyewall (due to a non-circular shape of the 
eyewall) that emerged from the boundary layer; these were classified with the class I/II 
trajectories by their starting points being below z = 1.5 km. The remainder of trajectories 
were then classified as class IV; thus these are trajectories with eyewall entry points above 
z = 1.5 km, radii exterior to the radius of maximum updraft, and starting locations above 
z = 1.5 km. 

3.4.2 Environmental entrainment 

Of the 400 trajectories, 213 are of class I/II, 115 are of class III, and 72 are of class IV. 
Class I/II trajectories are somewhat evenly distributed in azimuth (dark dots; Fig. 16b). 
The distribution of class III (triangles) and class IV (light dots; Fig. 16b) on the other 
hand show distinct azimuthal biases, apparently related to the 12 m s _1 northwesterly shear 
across the system; class III entry points appear on the downshear and class IV entry points 
on the upshear azimuths of the storm. The locations where the class IV trajectories have 
their maximum radius can be used (X’s in Fig. 16c, d) to study the means of interaction 
from the environment. These maximum radius points generally position to the left of the 
shear vector. Many of these trajectories quickly make their way to the eyewall, within a 
quarter rotation about the center; while others, whose maximum radius occurs up to five full 
hours prior to the seed time (the maximum possible by this calculation), orbit more than 
two times before encountering the eyewall. Several of those trajectories that make the quick 
approach to the eyewall appear to descend in approach along a slanted path, while the rest 
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appear to maintain a constant height in approach to the eyewall. 

3.4.3 Ventilation illustrated 

Ventilation is thought to weaken the hurricane, largely because of the transport of properties 
from the environment that are detrimental to maintenance of the hurricane. Our approach 
here is to investigate thermodynamic alterations of the eyewall by environmental air. The 
mid-troposphere in the tropical environment typically shows a mid-level 6 e minimum (Holton 
1992) ; in 0 e terms, the hurricane eye/eyewall are the warmest areas outside of the strato- 
sphere. Any mixing to the eyewall will reduce eyewall 6 e on average. A reduced 6 e has 
been hypothesized either to reduce the effectiveness of the hurricane heat engine (Riehl and 
Malkus 1961, Emanuel et. al 2004) and/or to result in a cooler warm core temperature 
anamoly aloft that supports the surface pressure signature of the hurricane (Gray 1968, 
Knaff et al. 2004) . 

The simplest measure of ventilation of eyewall air with environmental air is to note that 72 
of 400 trajectories representing 18% of the mass (and number, coincidentally) rising through 
z = 5 km in the eyewall originated in the environment (light colored dots, Fig. 16). The large 
number of interacting trajectories would appear to be a consequence of analyzing a sheared 
hurricane simulation. While a direct comparison of mixing processes in a three-dimensional 
simulation and axisymmetry must be done cautiously, we performed a comparative axisym- 
metric simulation using the 4x run from PM03 (which is neccessarily without vertical wind 
shear; not shown) and found zero class IV trajectories. We plan, in future work, to extend 
the present analysis to idealized, three-dimensional hurricane simulations without vertical 
shear. 

An analysis of entrant class IV trajectories is summarized in Table 1. Here, the change 
in 6 e along the trajectory is measured using 

A 9 e = 0 e (t = furthest radius) — 9 e (t = seed point) (1) 

and the dataset is partitioned according to whether A 6 e is positive or negative. When 
A 9 e < 0, the trajectory warms on approach to the eyewall, which must come at the cost of 
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a reduction of 9 e from the eyewall and the surrounding region, either through warming or 
moistening of the newly arriving trajectory. On average, class IV trajectories show negative 
values ({A 9 e ) = —0.89 K; where {•) is used for an ensemble average), indicating a net 
warming of such parcels upon approach to and entrance within the eyewall by some diffusive 
or irreversible mixing. Despite this warming, these class IV trajectories remain cool relative 
to the axisymmetric mean (9' e < 0). Also shown in Table 1 are standard deviations a of A 9 e 
and 9' e . Eyewall 9 e is reduced by the introduction of class IV trajectories, either directly by 
stirring cold air into the vicinity of the warmer eyewall air that emerges from the boundary 
layer inflow and the eye, or indirectly by diffusive mixing. Since class IV represents about 
18% of mass at z = 5 km, the axisymmetric mean eyewall is cooled 

0 18 

A 0 e (eyewall) « ___ (-0.89 - 3.47) K (2) 

« -0.96 K. 

The large standard deviations and the unbalanced sample size of positive and negative A9 e ’s 
are indicative of a positively skewed distribution, where, while all class IV parcels stir in 
mass, a small number of events may perform the largest modification of the eyewall mean 
Be- 

The technique of this section was applied using a second seed time (not shown) , t = 19 h 
instead of 20 h used above. At this time, the axisymmetric eyewall was narrower, so the 
seed annulus was 38 < r < 48 km. This provides 300 trajectories for analysis, 81 of which 
(27% ) are identified as class IV, thus there are more ventilating air parcels at this time. On 
the whole, this set of trajectories differs from the set illustrated above by a net decrease of 
entropy on approach to the eyewall ({A 9 e ) = 0.56 K instead of -0.89 K) with less variance 
(cr(A#e) = 3.86 K instead of 4.72 K), but a cooler asymmetric entropy at the seed point 
((B' e ) — —4.34 K instead of -3.47 K), and this set of trajectories approaches the eyewall more 
abruptly (change in azimuth (AA) = 4.15 rad versus 8.06 rad). Following the approach of 
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(3) above, the axisymmetric mean eyewall can be found at this time to cool 

0 27 

A0 e (eyewall) « (+0.55 - 4.34) K (3) 

w -1.40 K. 

The degree of eyewall cooling is found to be somewhat more consistent with the greater 
number of class IV trajectories. 

4 Summary and conclusions 

This study builds on the numerical study of Hurricane Bonnie by Braun et al. (2005) by 
examining the interaction between the various air masses (eye, eyewall, inflow, environment) 
of the TC. Through analysis of parcel trajectories, we have been able to -demonstrate that 
the low-level eye air can be carried to the eyewall; that this low-level eye air is replenished 
by inflow under the eyewall; that the mid- to upper-level eye air can be exchanged with 
eyewall air in a two-way process; and that mid- to upper-level air exterior to the eyewall can 
ventilate the eyewall (in the sense of Riehl and Simpson 1958, Riehl and Malkus 1961, and 
Emanuel et. al 2004). 

The trajectories examined demonstrate that the superintensity mechanism as described 
by PM03 and Braun (2002) is active in the Bonnie simulation. The results suggest that the 
superintensity mechanism is still able to benefit the intensity of the storm while it is being 
sheared by the environment. Trajectories originating in the low level eye and possessing 
relatively high 9 e are stirred into the base of eyewall updrafts and introduce an additional 
source of heat to the eyewall. The low-level eye is replenished by boundary layer parcels that 
slip underneath the eyewall updrafts to linger for some time in the eye (typically 40-60 min), 
thereby acquiring an elevated 6 e through prolonged interaction with the ocean surface. The 
stirring at low levels between the eye and eyewall occurs predominately left of the mean deep- 
layer. environmental shear vector. Our results are consistent with the analysis of updrafts 
in this MM5 simulation conducted by Braun et al. (2005), and are also supported by the 
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observational studies of Molinari et al. (2004) and Coborsiero (2004) and the theoretical 
studies by Reasor et al. (2004) and Jones (1995, 2004). 

More than half of trajectories originating inside the eye in the middle and upper tropo- 
sphere are not being contained inside the eye in our 5 h trajectory calculation period. More 
eye-eyewall mixing occurs in the upper troposphere than does in the hurricane boundary 
layer and middle (z = 5.6 km) troposphere. For mature hurricanes in low vertical wind 
shear conditions, W98 hypothesized that eye containment of air parcels above the eye inver- 
sion would be quite strong. Our case is more consistent with the weakening cases of W98 
and Carsey and Willoughby (2004). Containment of air within the eye is compromised in 
our case by stirring of air from the eyewall into the eye, but very little air is stirred into the 
eye from the environment above the boundary layer. Although a comprehensive analysis of 
mass replacement to the eye at mid- to upper levels was not carried out, there is evidence 
of eyewall air being stirred inward into the eye at these levels (e.g. Fig. 5). 

Ventilation (Riehl and Malkus 1961; Emanuel et. al 2004) is illustrated in the Bonnie 
simulation, as trajectories track the migration of environmental air with low 6 e into the 
eyewall. The net impact is about a I K reduction of eyewall average 0 e based on analysis of 
air at 2 = 5 km in the eyewall for two instances in time. For this simulation of a vertically 
sheared hurricane, the intensity is plausibly much greater than E-MPI (simulated 55 vs. 
E-MPI 38.2 to 52.2 m s -1 ). The interaction of environmental air with low entropy and the 
eyewall is hypothesized to weaken the TC. 

Future work will consider several of the questions raised herein in more detail. Specifi- 
cally, how substantial is 1 K cooling in the eyewall in reducing the intensity of the TC? A 
continued study of sheared storms in an idealized setting is ongoing to examine the efficiency 
of ventilation for various magnitudes of shear. Finally, it remains to identify the principal 
asymmetric structures around the storms that contribute the most to its ventilation. 
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Appendices 


A1 Determination of storm center 


The technique for finding the storm center described below basically attempts to place the 
maximum radial gradient in surface pressure close to the same radius from the center of the 
storm at all azimuths. The alternative method of using the location of minimum surface 
pressure was found to track centers of mesovortices rather than the storm-scale center of 
pressure. The technique described below is designed to avoid any arbitrary selection of scale 
(the implicit scale is the spatial scale over which the signature of the gradient of surface 
pressure is distributed) and any iterative processes for solution. 

The method for finding the storm center is as follows. Consider each grid point of a select 
subdomain of an x-y grid of simulated surface pressure. Consider each of four rays radiating 
along the four cardinal directions from each point. Obviously, grid points near the edge of 
the defined grid will lack data in one or two directions. For this reason, consider our select 
subdomain to be the middle l/9th of the domain, i.e., 


N x T 2 N x 

~ ^ ^ n 


(4) 

(5) 


where I x and I y are the index number of a grid point in the x- and ^-directions, and N x and 
N y are the total number of grid points in the x- and ^-directions, respectively. Determine 
the minimum M of the number of grid points lying along each ray. Along each ray, a profile 
of surface pressure Pj can be found, where j £ {1,2,3, 4} is an index for each of the four 
directions. Determine a “score” for each grid point s xy by a least mean-squares difference 
from their mutual average 


s = 


M 


( 6 ) 


where 



(7) 
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is the averaged pressure profile. Select the point where s is a minimum as the storm center. 

This method works because we know the innermost grid of the numerical model tracks the 
movement of the central surface pressure minimum, thus we are assured that the center will 
be found within the middle l/9th of the domain. Also, the scale of the vortex is somewhat 
larger than the scale of the inner-mesh grid. When the vortex scale is much smaller (such 
as on the coarsest grid mesh of our model run) than the size of the grid, then the four 
profiles might easily miss the vortex entirely and the score will measure differences in the 
more-or-less flat environmental surface pressure field, providing a small score trivially. 

We may take M above to be a fixed constant sufficiently small to handle every possible 
profile, i.e. 


M = min 


^ y. 


( 8 ) 


but we found that change did not alter the resulting centerfind for our case. 


A2 Computation of trajectories 

A trajectory is the computation of the path of a hypothetical fluid particle through a specified 
flow field 

-?(*„*) (9) 

where f m is the position of the trajectory of index m as a function of time and V is the 
Eulerian velocity field as a function of space and time. The trajectories are computed using 
a second-order Runge-Kutta method (Press et al. 1992) where four-dimensional (space-time) 
linear interpolation is used to compute values of the wind field that generally fall between 
spatial/temporal grid points. This simple problem is complicated by the shifting of the inner 
grid as it moves with the storm center. By modeling practice, the shifting of the grid moves 
not by one grid increment frequently, but by increments of several grid points occassionally. 
To solve this problem, a four-dimensional universal grid large enough to contain the complete 
translation of the grid as it follows the storm is created, and all information is remapped 
to this grid. Naturally, much of this universal grid is “empty”, and no effort is made to 
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populate these empty areas with other information available to us. Any trajectory leaving 
the inner (moving) mesh of the MM5 simulation is simply identified and considered to have 
left the neax-storm environment. 

The grid stagger of the wind fields is eliminated by interpolating the u, v, and w wind 
fields to the 9 or thermodynamic grid points, which is where the universal grid is defined. 
The winds are translated from m s -1 to the units of (gridpoints)s -1 , taking into account the 
effects of latitude using the formulae 


^grid 

^grid 


u u 

Ax cos (<f>) |AAjr@ 
v _ . v 

Ay (A4 >)r© 


( 10 ) 

(H) 


where u srid and u grid are wind speeds in grid units, (j) and A and latitude and longitude, 
respectively, converted to radians and r® is the mean radius of the Earth. Calculation in 
terms of grid-spacing per unit time is a coding convenience that allows a general-purpose, 
generic, well-tested trajectory routine to be used without being required to manage details 
of spherical geometry. If the geographic pole were included, we can expect the trajectory 
code to fail there unless a more appropriate model grid system were used. 

Since the trajectories are computed in Earth-relative coordinates, storm motion is not 
subtracted from the winds, but the resulting trajectories can be renavigated to storm- 
centered coordinates. 
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Figure 1: A simplified schematic of a moderately sheared hurricane. The dark (light) gray 
lines show the approximate outer (inner) edge of the eyewall. The dotted line represents 
the inversion level in the eye. The four subsections of section 3 examine different mixing 
process within the hurricane using different trajectory strategies for the numerical simulation 
of Hurricane Bonnie by Braun et al. (2005). Section 3.1 examines the mixing from the low- 
level eye to the eyewall (trajectory class I) by counting the number of trajectories that enter 
the eyewall versus those that remain in the eye. Section 3.2 examines the boundary layer 
inflow from the environment (class II) by considering the properties of trajectories that spend 
some time in the eye versus those that are transported immediately into the eyewall. Section 
3.3 examines the persistence of air in the middle and upper tropospheric eye (class III) by 
counting the number of trajectories that stir into the eyewall versus those that remain in the 
eye. Section 3.4 examines the ability of middle and upper tropospheric environmental air to 
stir into the eye or eyewall (class IV). 

Figure 2: Time series of minimum surface pressure (thick line), maximum tangential wind 
(thin line), and maximum azimuthally averaged tangential wind (dotted line) at the lowest 
model level (z = 40 m) for the period within the Hurricane Bonnie (1998) simulation of 
Braun et al. (2005) used here. The period covered here corresponds to 0800-1400 UTC 23 
August 1998. 

Figure 3: Contours of radial wind at (a) z — 453 m, (b) 1.1 km, (c) 5.6 km, and (d) 9.9 
km at t = 15 h. Positive values are shaded at levels of 0, 2.5, 5, 10 m s -1 . Negative values 
are contoured at -15, -10, -5, and -2.5 m s _1 . Dotted line represents the radius of maximum 
mean tangential wind at each level. 

Figure 4: Azimuthal averages of (a) tangential wind, (b) vertical velocity, and (c) equivalent 



potential temperature at t = 15 h in the 2 km Hurricane Bonnie simulation by Braun et 
al. (2005). Gray shading denotes azimuthally averaged relative humidity greater than 90%. 
The dotted line is the radius of maximum v as a function of height. 

Figure 5: Seven sample trajectories of class I (Fig. 1) viewed in (a) the "horizontal x-y plane 
and (b) the r-z plane. Trajectories are seeded less than 15 km from the storm center and at 
a height of 453 m. Trajectory seed time is 15 h into the 2 km grid simulation and trajectories 
are calculated 5 h forward in time. Associated values of 9 e as a function of height for the 
trajectories are shown in (c). Arrows indicate the direction of motion of the trajectories. 
The 5 h time average of the radius of maximum azimuthally-averaged tangential wind as a 
function of height is plotted as a dashed line in (b) . Seed locations are marked by light gray 
dots. 

Figure 6: Asymmetric (a) 9 e and (b) asymmetric w as a function of height along the trajec- 
tories of class I (Fig. 1) shown in Fig. 5. 

Figure 7: (a) Locations where trajectories seeded inside the eye ( 9 e > 370 K at a height of 
2 = 453 m (class I of Fig. 1) are stirred out to the axisymmetric eyewall. The arrow in the 
upper right corner indicates the direction of the mean layer (850 - 200 mb) vertical shear 
vector. Eyewall encounters occuring at or below z — 1.5 km are denoted by dark circles, 
encounters above z — 1.5 km are light gray, (b) Histogram of the azimuths of the points 
shown in (a). Of the 551 total trajectories seeded inside the eye at this level, 313 trajectories 
(57 %) are stirred out and meet the eyewall criteria. Eye seed points extend out to a radius 
of 24 km. 

Figure 8: (a) Histogram of trajectories of class II (Fig. 1) which are transported into the 
eye, shown as a function of trajectory seed radius. Trajectory seed points for this analysis is 
every 5 km in the radial direction and seed levels of z = 40 m, 121 m, and 244 m. Bin size 
is 5 km and sample size is 75 trajectories in each bin (25 at each seed level). Percentages 



are shown for (b) seed level z = 244 m, (c) seed level z = 121 m, and (d) seed level z = 40 
m. A trajectory is considered to be transported into the eye if its trajectory 9 e > 370 K at 
z <2 km. 

Figure 9: Radius versus height plots for a select sample of class II (Fig. 1) trajectories 
which can be identified as (a) “eye reservoir” trajectories or (b) “standard eyewall ascent” 
trajectories. 6 e versus height for the same trajectories are shown in (c) and (d) , respectively. 
Sample sizes are 41 for “eye reservoir” and 78 for “standard eyewall ascent.” 

Figure 10: Maximum 0 e versus eye residence time for the eye reservoir trajectories shown in 
Fig. 9(a,c). These trajectories are of class II from Fig. 1. 

Figure 11; Difference in 9 e between z = 1.5 km and 500 m (A 9 e = 9 e [z = 1.5 km] — 
9 e [z = 0.5 km]) for the “eye reservoir” (triangles) and “standard eyewall ascent” (squares) 
trajectories shown in Fig. 9(a,c). The numbering of trajectories along the x-axis is arbitrary. 
These trajectories are of class II from Fig, 1. • 


Figure 12: Skew T - log p diagram of the sounding located inside the eye at the storm center, 

/ 

The time displayed is t = 17.5 h, or halfway through the 5_h forward trajectory calculation. 
Temperature is the solid line and dew point is the dashed line. 

Figure 13: (a) Locations where trajectories seeded inside the eye (ui < 0.25 m s -1 and 
TL < 90%) at a height of z = 5.6 km are stirred out to the axisymmetric eyewall. The arrow 
in the upper right corner indicates the direction of the mean layer (850 - 200 mb) vertical 
shear vector, (b) Histogram of the azimuths of the points shown in (a). Of the 1497 total 
trajectories seeded inside the eye at this level, 878 trajectories (59 %) are stirred out and 
meet the eyewall criteria. These trajectories are of class III from Fig. 1. Eye seed points 
extend out to a radius of 40 km. 



Figure 14: (a) Locations where trajectories seeded inside the eye (w < 0.25 m s -1 and 
H < 90%) at a height of 2 = 9.8 km are stirred out to the axisymmetric eyewall. The arrow 
in the upper right corner indicates the direction of the mean layer (850 - 200 mb) vertical 
shear vector, (b) Histogram of the azimuths of the points shown in (a). Of the 1601 total 
trajectories seeded inside the eye at this level, 1210 trajectories (76 %) are stirred out and 
meet the eyewall criteria. These trajectories are of class III from Fig. 1. Eye seed points 
extend out to a radius of 42 km. 

Figure 15: Percentage of eye trajectories of class III seeded at each radius and at z = 453 m 
(solid line), 1.1 km (dotted line), and 5.6 km (dashed line) which satisfy the criteria for being 
stirred out to the axisymmetric mean eyewall (w > 0.25 m s _1 and H > 90%). Percentages 
are plotted as a function of trajectory seed radius. Sample sizes are 551 (453 m) , 903 (1.1 
km), and 1497 (5.6 km) trajectories. Eadial locations are binned to the nearest 1 km. 

Figure 16: (a) Locations where trajectories encounter the axisymmetric eyewall (w > 
0.25 m s _1 ) plotted in the (r, z)-plane. Dark circles are trajectories determined to be class 
I or II (BL inflow or low-level eye interaction; see Fig. 1); triangles are class III (mid- to 
upper- level eye interaction); and light circles are class IV (mid- to upper-level environmental 
interaction). The same is shown in plan- view in .(b). For each class IV trajectory, the eyewall 
entry point is again plotted with light circles and the point of furthest radius is shown with 
X’s in (c) an (r, 2 )-plane and (d) as a plan-view plot) with trajectories overlaid. The dashed 
line is the radius of maximum updraft. 
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Figure 1: A simplified schematic of a moderately sheared hurricane. The dark (light) gray 
lines show the approximate outer (inner) edge of the eyewall. The dotted line represents 
the inversion level in the eye. The four subsections of section 3 examine different mixing 
process within the hurricane using different trajectory strategies for the numerical simulation 
of Hurricane Bonnie by Braun et al. (2005). Section 3.1 examines the mixing from the low- 
level eye to the eyewall (trajectory class I) by counting the number of trajectories that enter 
the eyewall versus those that remain in the eye. Section 3.2 examines the boundary layer 
inflow from the environment (class II) by considering the properties of trajectories that spend 
some time in the eye versus those that are transported immediately into the eyewall. Section 
3.3 examines the persistence of air in the middle and upper tropospheric eye (class III) by 
counting the number of trajectories that stir into the eyewall versus those that remain in the 
eye. Section 3.4 examines the ability of middle and upper tropospheric environmental air to 
stir into the eye or eyewall (class IV). 



Figure 2: Time series of minimum surface pressure (thick line), maximum tangential wind 
(thin line), and maximum azimuthally averaged tangential wind (dotted line) at the lowest 
model level (z .== 40 m) for the period within the Hurricane Bonnie (1998) simulation of 
Braun et al. (2005) used here. The period covered here corresponds to 0800-1400 UTC 23 
August 1998. 
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Figure 3: Contours of radial wind at (a) 2 = 453 m, (b) 1.1 km, (c) 5.6 km, and (d) 9.9 
km at t = 15 h. Positive values are shaded at levels of 0, 2.5, 5, 10 m s -1 . Negative values 
are contoured at -15, -10, -5, and -2.5 m s -1 . Dotted line represents the radius of maximum 
mean tangential wind at each level. 





Figure 4: Azimuthal averages of (a) tangential wind, (b) vertical velocity, and (c) equivalent 
potential temperature at t = 15 li in the 2 km Hurricane Bonnie simulation by Braun et 
al. (2005). Gray shading denotes azimuthally averaged relative humidity greater than 90%. 
The dotted line is the radius of maximum v as a function of height. 
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Fig. 5: Seven sample trajectories of class I (Fig. 1) viewed in (a) the horizontal x-y plane 
and (b) the r-z plane. Trajectories are seeded less than 15 km from the storm center and at 
a height of 453 m. Trajectory seed time is 15 h into the 2 km grid simulation and trajectories 
are calculated 5 h forward in time. Associated values of 9 e as a function of height for the 
trajectories are shown in (c). Arrows indicate the direction of motion of the trajectories. 
The 5 h time average of the radius of maximum azimuthally-averaged tangential wind as a 
function of height is plotted as a dashed line in (b) . Seed locations are marked by light gray 
dots. 
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Figure 7: (a) Locations where trajectories seeded inside the eye ( 9 e > 370 K at a height of 
2 = 453 m (class I of Fig. 1) are stirred out to the axisymmetric eyewall. The arrow in the 
upper right corner indicates the direction of the mean layer (850 - 200 mb) vertical shear 
vector. Eyewall encounters occuring at or below z = 1.5 km are denoted by dark circles, 
encounters above 2 = 1.5 km are light gray, (b) Histogram of the azimuths of the points 
shown in (a). Of the 551 total trajectories seeded inside the eye at this level, 313 trajectories 
(57%) are stirred out and meet the eyewall criteria. Eye seed points extend out to a radius 
of 24 km. 





Figure 8: (a) Histogram of trajectories of class II (Fig. 1) which are transported into the 
eye, shown as a function of trajectory, seed radius. Trajectory seed points for this analysis is 
every 5 km in the radial direction and seed levels of z = 40 m, 121 m, and 244 m. Bin size 
is 5 km and sample size is 75 trajectories in each bin (25 at each seed level). Percentages 
are shown for (b) seed level z — 244 m, (c) seed level z = 121 m, and (d) seed level z = 40 
m. A trajectory is considered to be transported into the eye if its trajectory 9 e > 370 K at 
2 < 2 km. 
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Figure 9: Radius versus height plots for a select sample of class II (Fig. 1) trajectories 
which can be identified as (a) “eye reservoir” trajectories or (b) “standard eyewall ascent” 
trajectories. 0 e versus height for the same trajectories are shown in (c) and (d), respectively. 
Sample sizes are 41 for “eye reservoir” and 78 for “standard eyewall ascent.” 
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Figure 10: Maximum 6 e versus eye residence time for the eye reservoir trajectories shown in 
Fig. 9(a,c). These trajectories are of class II from Fig. 1. 
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Figure 11: Difference in 6 e between z = 1.5 km and 500 m (A 9 e = 6 e [z = 1.5 km] — 
6 e [z = 0.5 km]) for the “eye reservoir” (triangles) and “standard eyewall ascent” (squares) 
trajectories shown in Fig. 9(a,c). The numbering of trajectories along the x-axis is arbitrary. 
These trajectories are of class II from Fig. 1. 
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Figure 12: Skew T - log p diagram of the sounding located inside the eye at the storm center. 
The time displayed is t = 17.5 h, or halfway thrpugh the 5 h forward trajectory calculation. 
Temperature is the solid line and dew point is the dashed line. 
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Figure 13: (a) Locations where trajectories seeded inside the eye ( w < 0.25 m s _1 and 
H < 90%) at a height of 2 = 5.6 km are stirred out to the axisymmetric eyewall. The arrow 
in the upper right corner indicates the direction of the mean layer (850 - 200 mb) vertical 
shear vector, (b) Histogram of the azimuths of the points shown in (a). Of the 1497 total 
trajectories seeded inside the eye at this level, 878 trajectories (59 %) are stirred out and 
meet the eyewall criteria. These trajectories are of class III from Fig. 1. Eye seed points 
extend out to a radius of 40 km. 
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Figure 14: (a) Locations where trajectories seeded inside the eye (w < 0.25 m s -1 and 
H < 90%) at a height of z = 9.8 km are stirred out to the axisymmetric eyewall. The arrow 
in the upper right corner indicates the direction of the mean layer (850 - 200 mb) vertical 
shear vector, (b) Histogram of the azimuths of the points shown in (a). Of the 1601 total 
trajectories seeded inside the eye at this level, 1210 trajectories (76 %) are stirred out and 
meet the eyewall criteria. These trajectories are of class III from Fig. 1. Eye seed points 
extend out to a radius of 42 km. 





Figure 15: Percentage of eye trajectories of class III seeded at each radius and at z = 453 m 
(solid line), 1.1 km (dotted line), and 5.6 km (dashed line) which satisfy the criteria for being 
stirred out to the axisymmetric mean eyewall (w > 0.25 m s -1 and TL > 90%). Percentages 
are plotted as a function of trajectory seed radius. Sample sizes are 551 (453 m), 903 (1.1 
km), and 1497 (5.6 km) trajectories. Radial locations are binned to the nearest 1 km. 
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Figure 16: (a) Locations where trajectories encounter the axisymmetrie eyewall (w > 
0.25 m s -1 ) plotted in the (r, z)-plane. Dark circles are trajectories determined to be class 
I or II (BL inflow or low-level eye interaction; see Fig. 1); triangles are class III (mid- to 
upper-level eye interaction); and light circles are class IV (mid- to upper-level environmental 
interaction). The same is shown in plan-view in (b). For each class IV trajectory, the eyewall 
entry point is again plotted with light circles and the point of furthest radius is shown with 
X’s in (c) an (r, V) -plane and (d) as a plan-view plot, with trajectories overlaid. The dashed 
line is the radius of maximum updraft. 






Popular Summary 


A Lagrangian trajectory view on transport and mixing processes between the eye, eyewall, and 
environment using a high-resolution simulation of Hurricane Bonnie (1998). 

By T. A. Cram, J. Persing, M. T. Montgomery, and S. A. Braun 
Submitted to the Journal of the Atmospheric Sciences 

Hurricane intensity is controlled by many factors, including the characteristics of a 
storm’s environment as well as processes internal to the storm. The theoretical maximum 
potential intensity that a hurricane can achieve is a function of the underlying sea surface 
temperature (SST) and the temperature of the exhaust outflow at high levels of the troposphere, 
with warmer SSTs and colder outflow temperatures leading to potentially stronger storms. The 
effects of environmental vertical wind shear on a hurricane are thought to transport 
environmental air into the storm, a process sometimes referred to as ventilation, and weaken its 
warm core and strong winds. Cold downdrafts and air-sea interactions associated with storm- 
induced cooling of ocean temperatures also tend to reduce the ability of a storm to reach its 
maximum potential intensity. A hurricane has a primary circulation, which is the strong swirling 
flow around the center, and a secondary circulation, which consists of a radially in-up-and-out 
flow pattern which serves to transport heat to the upper troposphere. This in-up-and-out 
circulation plays a major role in the intensification process and is central to the theory of 
maximum potential intensity. An interesting characteristic of hurricanes simulated in ideal 
conditions is that, if model resolution is fine enough to resolve mixing processes in the eyewall 
region, the simulated storms are capable of reaching intensities greater than the theoretical 
maximum value. This state is referred to as “super-intensity”. Previous studies have shown that 
mixing of air between the hurricane eye and the eyewall may act similar to a turbo-booster by 
injecting high-energy air in the eye into the eyewall to produce stronger updrafts and a stronger 
secondary circulation. 

In this study, we use a high-resolution simulation of a hurricane that was observed in 
1998 during the NASA Convection and Moisture Experiment (CAMEX-3) to study mixing 
processes in the storm and their impact on storm intensity. Using output from the model, the 
hurricane is seeded with thousands of virtual particles that move with the storm’s winds. These 
particles are followed around the storm and are used to characterize mixing between the eye and 
eyewall as well as between the eyewall and its near environment outside of the eyewall. The 
trajectory analysis shows that 

• Low-level eye air— possessing the highest thermodynamic energy in the storm— is mixed 
into the eyewall, provides an additional source of heat to the eyewall updrafts, and helps the 
storm to reach or slightly exceed its maximum potential intensity despite the negative effects of 
vertical shear and cold downdrafts. 

• The low-level eye is replenished by inflowing boundary layer air parcels that slip 
beneath the eyewall and linger for some time in the eye, thereby acquiring greater 
thermodynamic energy. 

• About 20% of the air at middle levels within the eyewall has its origins in the middle to 
upper level environment, which is characterized by lower thermodynamic energy. This lower 
energy air is introduced into the eyewall by the ventilation effects of the environmental wind 



shear and produces a significant cooling and/or drying effect, which helps to limit the storm’s 
intensity. 

In 1995, we entered an active phase of a multi-decadal hurricane cycle in the Atlantic 
Ocean, meaning that the United States will remain vulnerable to hurricanes for many years to 
come. While forecasts of hurricane tracks have improved over the past thirty years, forecasting 
of storm intensity remains a difficult problem due to a combination of inadequate observations 
and models as well as a poor understanding of the physical and dynamical processes involved. 
The present study is important because it sheds light on some of the key processes that influence 
intensity. 



